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Abstract. The foundations of weak turbulence theory is explored through its ap- 
plication to the (alpha) Fermi-Pasta-Ulam (FPU) model, a simple weakly nonlinear 
dispersive system. A direct application of the standard kinetic equations would 
miss interesting dynamics of the energy transfer process starting from a large-scale 
excitation. This failure is traced to an enforcement of the exact resonance condi- 
tion, whereas mathematically the resonance should be broadened due to the energy 
transfer happening on large but finite time scales. By allowing for the broadened 
resonance, a modified three-wave kinetic equation is derived for the FPU model. 
This kinetic equation produces some correct scaling predictions about the statistical 
dynamics of the FPU model, but does not model accurately the detailed evolution of 
the energy spectrum. The reason for the failure seems not to be one of the previously 
clarified reasons for breakdown in the weak turbulence theory. 

1. Introduction. Weak turbulence (WT) theory was developed in the 1960s to 
provide equations which quantitatively describe the transfer of energy among tur- 
bulent, weakly nonlinear, dispersive waves in fluids |^, |T[ The 
basic "kinetic equations" produced by WT theory have been applied to analyze en- 
ergy transfer include internal and surface waves with small aspect ratios in the at- 
mosphere and ocean , semiconductor lasers , and plasma turbulence . 
The domain of validity of the governing equations, however, remains to be fully 
explored. 

Zakharov and Kadomtsev obtained the kinetic equations through ar- 
guments (such as random phase approximation) motivated by physical intuition, 
whereas various authors in the USA and Western Europe ^ ^ derived the 
same kinetic equaitons through a multiple scales expansion. While the latter ap- 
proach is systematic, it has not been rigorously formulated, and several breakdowns 
of its validity have been noted in recent years. Newell and coworkers Q ||, |2l| have 
noted that the basic asymptotic expansion parameter is usually not uniformly small 
over the full range of scales in a weakly turbulent system, so that "intermittent" 
effects not directly predictable by WT theory will manifest themselves at either suf- 
ficiently large or small scales. Cai, Majda, McLaughlin, and Tabak |^, |l9j have 
scrutinized the validity of the kinetic equations of WT theory and their underlying 
assumptions through their application on a simple one-dimensional model amenable 
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to comparison with direct numerical simulation. One key finding of their work is 
that under certain circumstances, a nonlocal coupling between the large-scale forc- 
ing and much smaller "inertial range" scales can arise which invalidates the tacit 
assumption that the energy transfer in this inertial range of scales is unaffected by 
the presence of forcing and/or damping. 

To further these efforts at understanding the domain of validity of WT theory, we 
study its application to a simple Hamiltonian model without driving or damping. 
The Fermi-Pasta-Ulam (FPU) model of a one-dimensional collection of particles 
coupled by weakly anharmonic springs serves as a natural test example which osten- 
sibly should be amenable to the WT approach, provided the number TV of particles 
in the chain is large enough. Indeed, the FPU model has dispersive waves which 
are coupled through a weak nonlinearity, and it can be initialized with a disordered 
wave profile. Moreover, its Hamiltonian nature excludes the problem related to the 
nonlocal force coupling. We consider here the a version of the FPU model (defined 
in Section ^ which has a cubic nonlinearity in the Hamiltonian and therefore direct 
three-wave coupling (which is usually simpler to analyze than four-wave couplings 
as studied in |7[ |l^). Nonetheless, we still find the standard kinetic equations of 
weak turbulence insufficient to predict the dynamics of the energy spectrum in the 
FPU model when it is initialized with a large-scale excitation. 

As we discuss in Section ^ after some notational preparation in Section |^, there 
are no nontrivial exact triad (three- wave) resonances. A literal application of the 
standard WT theory would suggest then that the energy transfer between waves 
must be governed by four-wave interactions. A scrutiny of the numerical simula- 
tions of the Qf-FPU dynamics however, indicates the presence of at least two 
distinct time scales on which significant energy transfer occurs. The first phase 
only involves energy transfer and equilibration among modes of sufficiently small 
wavenumber. In the final phase, occurring on a much longer time scale, energy is 
redistributed among all modes. These observations can be reconciled with the reso- 
nance properties of the dispersion relations by recognizing that the low wavenumber 
modes are in near resonance because the waves become approximately acoustic at 
low wavenumber. Consequently, we can interpret the first phase of energy transfer 
among low wavenumber modes as being mediated by triad interactions and the later 
phase of relaxation to global energy equipartition as requiring quartet interactions. 

We attempt here to adapt weak turbulence theory to describe the energy transfer 
among nearly resonant triads during the earlier phase of evolution of the energy 
spectrum in the a-FPU model. Revisiting the multiple scales derivation of the 
three- wave kinetic equation in Section ^, we identify a step in which a large but 
finite time is replaced by an infinite time limit. It is this simplification which enforces 
the exact resonance condition present in the standard kinetic equations. In reality, 
the resonance should be broadened simply because the three-wave kinetic equation 
seeks to describe the statistical dynamics of a turbulent system on at a large but 
finite time. Over a long but finite time T, frequency mismatches of order o(T~^) 
cannot be detected. By properly retaining the effects of this "finite-time resonance 
broadening" and employing a systematic renormalization group technique, we can 
compute the range of wavenumbers which will exchange significant energy through 
triad interactions and the time scale on which this energy exchange should occur. 
These scaling predictions are well confirmed by the numerical results, and agree 
with earlier partial predictions |ll|, psj . 
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Our calculations also produce a modified kinetic equation ( |12| ) with finite-time 
resonance broadening. Unfortunately, the details of this kinetic equation are not in 
good agreement with the numerically observed energy spectrum evolution. More- 
over, the modified kinetic equation generates negative energies. The problem seems 
to be the omission of higher order secular terms in the multiple scales expansion 
which play an important role when the waves are weakly dispersive [ p^ , but does 
not appear to fall neatly into an established situation where the WT theory is known 
to fail. In future work, we will attempt to clarify the nature of the breakdown and 
assess whether a further modified WT approach might succeed. 



2. Fermi-Pasta-Ulam Model. The Fermi-Pasta-Ulam (FPU) model is a model 
for a one-dimensional collection of particles with massless, weakly anharmonic (non- 
linear) springs connecting them to each other. Letting {qj}jLi and {pj}jLi denote 
the position and momentum coordinates of an A'^-particle chain, we can define the (a 
version) FPU model Hamiltonian, which we write immediately in nondimensional 
form: 

N N 

pj 

[l]-! - 2gj + qj+i) (1 + e(gj-i - qj+i)) 

The coordinates are normalized so that the total Hamiltonian energy H = N, 
implying a typical 0{1) value of each qj and pj. The small parameter e <C 1 then 
measures the weakness of the nonlinearity. 

We now make a conformal transformation to action-angle variables to facilitate 
the application of WT theory. We describe the transformation by the following two 
steps: 




where uJk = 2|sin(^)| is the dispersion relation describing the frequency of a 
discrete Fourier mode k. The virtue of the new (complex- valued) canonical variables 
ak is that, in the absence of nonlinearity, they rotate with a single frequency ( 
ak (X e"*"^*"*). To connect with the standard notation of WT theory and keep 
subsequent expressions as compact as possible, we define the variables A'^ = a^, 
and = alj.. The superscript refers to the branch of the dispersion relation: 



n = 

dq^ 

dt 
dt 
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In these new variables, the Hamihonian reads 



N 



^ = ^ E E ^r;i:^^i^r.4.+fc..o (3) 

fel,fc2 = -f + 1 «1.S2 = ± 

N_ 
2 

I ^ \ ^ \ ^ n_lSi,S2,S3 ASi AS2 AS3 X 

fcl,fe2,fc3 = --f +1 Sl,S2,S3 = ± 

Wfeia = \v^w;5s,.-s2, (4) 



^si,^2,S3 _ 51 52 S3 Sgn (fci fca fc3 ) ^U}k^ ^fc^ ^fc3 

The symbol Sij denotes the Kronecker delta function. The equations of motion can 
then be written [p2[ 

A AS 1 2 

-^-^su;al = E E ^ra^".^r2'^fc.'=.+^- (6) 

^ssi^, ^ i5i52 sgn (fcfcifcz) yg^fcC^fciC^fc^ 
kkik2 ' 

To keep focus, we will restrict the initial data to have statistical symmetry under 
parity inversion and to have Gaussian statistics. Then it suffices to describe the 
energy transfer among different scales through the statistical quantities |^ : 

nk^l^iH?) = ^(A+A:,), (8) 

Ek = uJkTik + LU^kn{-k) = 2uJknk. 

The main object of our interest will be Ek, the harmonic energy associated to the 
wavenumber ±fc, or in other words the mean contribution to the quadratic part 
of the Hamiltonian in (|^) from wavenumbers ±fc. The mean-square amplitude Uk 
of the action variables is sometimes referred to as the "occupation number" of the 
mode k. 



3. Notational Convention. The weak turbulence theory which we will apply to 
analyze the FPU model is usually formulated in terms of spatially homogenous 
systems of infinite extent, so that the Fourier wavenumbers are continuously dis- 
tributed. This would be appropriate for a particular well-defined N oo limit 
of the FPU model, in which the physical lattice spacing remains order unity and 
the system becomes infinitely extended in both directions in physical space. We 
will always consider N to be finite but large, and introduce notation which makes 
contact with that usually used in the weak turbulence literature j22|. 

We define the normalized wavenumber ^ = k/N which takes values in the unit 
circle Si (or equivalently, the space of real numbers M with all numbers differing 
by an integer declared equivalent) To distinguish between Fourier coefficients 
defined with respect to the integer wavenumbers k and the new wavenumber variable 
^, we write Ofc for the former and a(^) for the latter. The associated frequency is 
defined through the continuous dispersion relation lu{S,) = 2 |sin7r^|. 
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For large but finite N, the correlation functions of the Fourier coefficients are 
proportional to a discretized delta-function, in the sense that 

(a(o«*(n) = n{o~s^^^H^-a 

where n{^) has a finite, nontrivial iV — > oo limit. 

To further emphasize the connection between the finite but large N and infinitely 
long systems, we introduce the discretized integral notation: 

which should just be thought of as a Riemann integral approximation to the integral 
defined with respect to a mesh situated at integer multiples of 1/M. 

4. Modified Kinetic Equation. The standard weak turbulence theory |2^, ^ 
would predict that three-wave interactions would be incapable of significant en- 
ergy transfer. The reason is that the only triads of waves which satisfy the exact 
resonance conditions 

? = a + 6, c^(o = ±c^(a)±^(6) 

are those involving the zero mode = or = or ^2 = 0), which is uncoupled 
to the other degrees of freedom due to the translational symmetry of the system 
(-^5516 = if ^^1^2 — 0)- The absence of three- wave resonances is here not due to 
the finite size of the system but simply due to the concavity of the dispersion 
relation Q. 

The reason the standard kinetic equations of weak turbulence theory fail to de- 
scribe any energy transfer through three- wave interactions is that these interactions 
are effective only over a thin band of wavenumbers ^ with width vanishing with the 
nonlinearity parameter e. The derivation of the standard weak turbulence kinetic 
equations assumes that all quantities, including the wavenumbers ^ of interest and 
the energy density E{^) in each mode, are order unity in magnitude (asymptotically 
independent of e). The fact that the wavenumbers of interest during the three- wave 
dynamics are small, and the energy density residing in these wavenumbers large, 
requires a reworking of the asymptotic derivation of the kinetic equation. 

To this end, a simple renormalization group approach is useful . We rescale 

the wavenumber and amplitude of the energy spectrum by some gauge functions 
/3(e) and 7(e) depending on s: 

e = e//3(e), t = t/jie). (9) 

The rescaling factors f3 and 7 are chosen so that ^ and t are order unity for the 
three-wave dynamics of interest. We shall leave them unspecified for now, other 
than insisting that /? J, and 7 — > cxd as e ^ 0, and show how their scaling with e 
can be determined by self-consistency. More specifically, we define a rescaled energy 
spectrum i?(^)(^, t), obtained by by applying the transformation (^ to the original 
energy spectrum: 

EiU)=P{e)-'E(^\am,tMe)) (10) 
The rescaling of the amplitude results from the normalization condition that the 
total quadratic energy is N, which implies that if E{^, t) is concentrated over a 
thin band of wavenumbers |^| < /^(e), then the average amplitude of the energy 
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spectrum over these modes must be oid{(3{e)~-^). The rescaled energy spectrum 
-E*-'^-' t) is supposed to be a truly order unity function, with order unity amplitude 
and order unity variations with respect to ^ and t during the three-wave dynamical 
phase for all sufficiently small e. Mathematically, we demand that the renormalized 
energy spectrum 

be a nontrivial function with order unity amplitude and derivatives. In renormal- 
ization group terminology, we want E{^,t) to be a nontrivial fixed point under the 
transformation induced by (0). The gauge functions /3(e) and 7(e) are determined 
by this property. Indeed, as we shall show in Section 

/3 = ei/2, 7 = £-3/2. (11) 

and the renormalized energy spectrum satisifes the following modified kinetic equa- 
tion: 

X [p(^"i,t)^(^"2,t) - ^iEiU)m2,t) - ^2m,t)mi,t)] da 

(13) 

We comment now on these predictions. 

4.1. Scaling Predictions. The self-consistency calculation based on the weak tur- 
bulence expansion predicts that the band of wavenumbers participating in the ear- 
liest (three-mode) phase of energy transfer has thickness Afc ^ (3 — e^/^ and asso- 
ciated dynamical time scale i ~ 7 = e"'^/^, provided the initial data is distributed 
over a band of low wavenumbers which is comparably (e^^^) thick. In these 
scaling predictions are confirmed through direct numerical simulations. Analytical 
explanations for the e^/^ width for the band of wavenumbers which equilibrate first 
have been previously offered on the basis of resonance overlap calculations [p5| and 
related considerations |ll|. One can also construct an argument for the e^'^/^ 
dynamical time scale from cruder arguments. Our intent, however, was to provide a 
more systematic argument, since resonance overlap arguments can sometimes give 
misleading results |15|] . As we discuss below, however, we cannot yet claim to have 
fulfilled this goal. 

4.2. Properties of Modified Kinetic Equation. The modified kinetic equation 
T2|) has the following good properties: 

The modified three-wave kinetic equation ( p^ conserves the harmonic energy 
over the three-wave time scale: 

More precisely, direct calculations on the modified kinetic equation ( p^ show 
that the three- wave interactions conserve the harmonic energy locally (triad by 
triad) for the band of wavenumbers ^ < e^/"^ ^ while the energy in wavenumbers 
^ ^ e^/2 is predicted not to change significantly on the time scale t < e"'^/^. 
The modified kinetic equation accounts for a finite time resonance broadening 
in that the usual Dirac delta function S{Auj) involving the frequency differ- 
ence ( |l6|) is instead replaced in (^2|) by a renormalized form of the function 
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'^'"aL"*^ • While the latter does indeed converge to the Dirac delta function 
in the t oo limit, its distinction from this strict limit is important for pre- 
dicting nontrivial dynamics for large but finite t. Some previous work ^ 
offered ad hoc modifications to the weak turbulence theory to take into ac- 
count such broadening of the resonances, but we have endeavored to derive the 
quantitative form of the finite-time resonance broadening through systematic 
means. 



4.3. Deficiency in Modified Kinetic Equation. The modified kinetic equation, 
however, has a serious deficiency in that it does not preserve nonnegativity of the 
energy spectrum, because the broadened resonance kernel sin(7r'^ffi^2i)/(ffiC2) is 
not sign-definite. This indicates that, unlike standard weak turbulence theory, the 
removal of the leading order secular term (along with frequency renormalization) 
does not simultaneously remove all higher order secularities active at the first non- 
linear time scale, t ^ ^ = £~^/^, in the a-FPU model. A similar issue arose in the 
problem of acoustic turbulence pq | . 

Newell and coworkers ^ |21| have recently pointed out that the assumptions 
underlying the WT theory will often fail at either large or small scales, because the 
nonlinear time scale tNL associated with such scales may fail to be long compared 
with the corresponding linear time scale tj^. Because the development of partial 
equipartition is concentrated on small wavenumbers k < P ^ e^^'^, one might think 
that the deficiencies noted for the kinetic equation derived above arise from this 
kind of breakdown in the asymptotic expansion. We now show, however, that the 
simplest such explanation is insufficient. 

The linear time scale th associated with a wavenumber ^ is just the inverse 
frequency while the associated nonlinear time scale iNL ^^^jTf^- Fi'om 

studying either the primitive second order WT expansion ( [1^ ) or the renormalized 
kinetic equation, we see that for the band of active wavenumbers ^ < e^/^, the 
nonlinear time scale depends on f and e as <nl Therefore, for ^ < e^/^, 

the ratio ^nl/^l ~ e"^ is uniformly large, so the defects in the kinetic equation 
cannot be attributed to this criterion elaborated in jj, |9[ ^ . 

Another situation in which weak turbulence theory will fail is in finite systems 
where the wavenumber spacing is not small compared with the width of the 
frequency resonance, which is here proportional to /3 ~ e^^^ Therefore, WT 

theory cannot be expected to apply to the a-FPU model when Ne^^^ < 1, but this 
can be avoided, at least in principle, by choosing the size of the system N to be 
large enough. 

At this time, we are endeavoring to specify concretely the source of the deficiency 
of the modified kinetic equation. 



5. Derivation of Modified Kinetic Equation. Our starting point will be the 
following equation, which results from a multiple-scales expansion for the moments 



8 



PETER R. KRAMER, JOSEPH A. BIELLO, AND YURI V. LVOV 



of the mode amplitudes a ^ 



MN) MN) 

= -4e'n / / /C(*)(c^(e)-sic.(ei)-.S2c^(6)) 

x^«W(C - a - e2)e^^-"(«^+^^"(«^^+^^"(«^»* d^ dC2 



dt 



+ ... 



(14) 



Here we have used the symbol 3? to denote the real part of the ensuing expression, 
and defined the function 



1 



The ellipses denote terms which are omitted because they are formally of higher 
order in e. We will revisit the effects of these neglected terms in future work, because 
terms which are insignificant in standard weak turbulence theory must be carefully 
examined to see if they remain insignificant in our modified derivation. 

We convert Eq. (^) into an equation for the energy spectrum through the relation 
E{^) = 2ti>(^)n(^) and the explicit formula (|^) for the coupling coefficient. Then 
upon substituting the rescahng (p^) into we obtain 



dt 



.„..„=+ -'/J-^Si 



Si,S2=± 



,E('Hh,t)E(^Htt) 



_ f — SljS2,- 

'^-/35i,/3f2,-/3C' 



-Lj^:jilE^'^ilt)E('H^i,t) 



dCi d6 



where 



and 



(15) 



(16) 



rSSiS2 



I ^iO rssiS2 ^ isiS2 sgn(fcfcifc2)cj(0 



We now seek particular relations between 7, /3, and e so that the right hand side of 
(^5|) is consistent with the assumed order unity magnitude and scale of variation of 
the left hand side in the limit e ^ 0. To this end, we recall that /3(e) | 0, 7(e) 00, 
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and observe that 



Aw 



2nm\-si\ii 
^/3(IC|-Sil6|-^*2|C2 



-icKsiiai 



- 52 16 



52i6n 



el 



el 



siS2iaii6i) 



Noting that ^ = 6 + 6 from the Kronecker delta function in (|15|), the appropriate 
asymptotics for Aw can be decided based on the si, §2, and the signs of 6 61 
and 6- If the signs ahgn so that |6 + 61 = si|6l + S2|6|j then Aw ^ ord(/3'^); 
otherwise Aw '--^ ord(/3). Moreover, 




i'r3siS2l«ll«lll«2l 



2^J(|«|-Sl|«l|-S2|«2|) 



0(/32) 



if 16+61 

otherwise 



S1I6I + S2I6I 



Therefore, the dominant contribution will come from the wavenumbers in the do- 
main such that 16 — 16+61 = S1I6I + 52161- We will choose the scaHng functions 
so that Aw < 7^^ for order unity values of 6 61 ^^'^ 6 in this dominant region, 
so we can simplify further: 

e-'^'^"*/C®(7Aw) = 7^^/3-3 

when 16 



giTr^7/3^siS2|5||5i||52|t _ I 

*^3siS2l6l6ll6l 

+ 61 =Sl|6l+ 52161- 



Discarding subdominant contributions based on the above considerations, we 
achieve after some analytical manipulations: 



lim 



dt 



lim 



£7 



-0 TT(3{e) 



"f'^^) /■«''^) sin(^V!|6M 

66 



X p(6, f)^(6, - 6^(6 i)^(6, - 6^(6 i)^(6, 
x5»(/3^)(e-6-6)d6d6 

The equating of the limit of the integral on the right hand side by the integral of 
the limiting integrand can be justified through uniform bounds on finite intervals 
and the relative negligibility of the contribution from large 6 ^nd 6- The latter 
fact can be checked by direct analysis or more intuitively understood from the 
observation that the wavenumbers with ord(l) magnitude ^ (and therefore large 
values of f = exhibited no activity on the three-mode interaction time scale. 

The gauge functions /3(e) and 7(e) were introduced so that E'^'^^ (6 1) would have 
a nontrivial dynamics with order unity magnitude in the e — *■ limit. We see then 
that a nontrivial, finite dynamics results precisely if 0^^ and e'^jl3~^ have finite, 
nonzero limits. We thereby arrive at the unique (up to irrelevant multiplicative 
constants) choice of gauge functions ( |ll|) and modified kinetic equation (|l^) for the 
renormalized energy spectrum. 
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